The moment sum rule and its consequences for ferromagnetism 

in the Hubbard model 

M. Potthoff, T. Herrmann, T. Wegner and W. Nolting 
Lehrstuhl Festkorpertheorie, Institut fur Physik, Humboldt- Universitat zu Berlin, Germany 



The sum rule for the moments of the spectral density is discussed for the single-band Hubbard 
model. It is shown that respecting the sum rule up to the order m = 3 is conceptually impor- 
tant for a qualitatively correct description of the quasi-particle band structure in the strong- 
correlation regime. Different analytical approximations for the self-energy are analyzed with 
respect to their compatibility with the moment sum rule. To estimate the practical usefulness 
of the sum rule, correlation functions and dynamical quantities are determined. The results 
obtained within the various approximation schemes of different complexity are compared 
with each other and also with essentially exact results available for infinite-dimensional lat- 
tices. It turns out that the m = 3 moment is rather unimportant for the paramagnetic phase 
on the hyper-cubic lattice. Contrary, it decisively influences the magnetic phase boundary 
as well as the critical temperature for the ferromagnetic phase on an fec-type lattice. 
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I. INTRODUCTION 

The Hubbard model [|, ||, ||] plays a central role in 
the attempts to understand the physics of correlated 
electrons on a lattice. It is surely oversimplified when it 
comes to a description of real materials such as the 3d 
transition metals and their oxides, for example. Never- 
theless, it can provide deep insight into the fundamental 
mechanisms that are responsible for various prominent 
correlation effects. 

The Hubbard model sets up a notoriously difficult 
problem which now as before is not accessible to an ex- 
act solution in general. An exception is given for the 
one-dimensional case (d = 1) where a yet incomplete 
but very detailed understanding of the model proper- 
ties in the whole parameter range has been achieved 
@' H' H' 0' H ■ An important simplification of the model 
is also given in the opposite limit of high spatial di- 
mensions d = oo |9| [l0[ [ll], |l2| which was recognized 
as a well-defined and non-trivial limiting case by Met- 
zner and VoUhardt B. Here the electronic self-energy 
is k independent or site-diagonal and thereby the 
model becomes equivalent to an effective impurity prob- 
lem |Q . An exact solution is possible by self-consistent 
mapping jn| [l6| onto the single-impurity Anderson 



model (SIAM) fl7||, for example, followed by the numer- 
ical treatment [18] within the Quantum-Monte-Carlo 
(QMC) approach |l6|, (jJI, §f| In dimensions d = 3 
or d = 2, however, one still has to resort to approximate 
treatments. 

Valuable information that helps to judge of the reli- 
ability of a particular method is provided by exact iden- 
tities, sum rules, limiting cases etc. In any dimension 
d such rigorous results impose strong necessary condi- 
tions for the inevitable approximations. The purpose 
of the present paper is to focus on a particular sum 
rule: Exact expressions can be obtained for the (low- 
order) moments J E m A\ La {E) dE of the spectral den- 
sity Ak a (E). We will argue that the first four moments 
(m = — 3) yield valuable information on the quasi- 
particle band structure; they are especially important in 
the strong-coupling regime and also decisively influence 
the possibility and characteristics of spontaneous mag- 
netic order. The moment sum rule has been considered 
not only within the context of the standard single-band 
Hubbard model in d = oo pi pi §, §, in d = 3 
|,§§|,| l( | = 2 |lM H and d = 1 |3 , but 
also for the negative-C/ case J34J and for reduced trans- 
lational symmetry [ff, |§, for the SIAM fl, for the 
t-J Q and for localized spin models |3§, 4C, and 
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may thus be of general interest for the construction of 
analytical approaches. 

There are several questions related to the moment 
sum rule which are not yet finally clarified. Firstly, we 
have to ask how to check to which order the sum rule 
is fulfilled for a particular (approximate) method. Sec- 
ondly, a kind of recipe is required that shows up how a 
given method can be modified to respect the sum rule up 
to a certain desired order. Thirdly, it should be worked 
out what is the actual conceptual improvement that is 
thereby achieved. Finally and most important, we need 
to know what can be achieved in practice, i. e. whether 
the sum rule can help to come "closer" to the exact so- 
lution. 

Especially the last question requires to compare with 
exact results. As has been mentioned above, these are 
available for the Hubbard model in infinite dimensions. 
In the recent years there has been extensive work for 
d = co lattices concerning the Fermi-liquid phase at 
and off half-filling |2I[| , the metal-insulator (Mott) tran- 

transport properties 

II, m 6 H 

9|. These QMC 



sition 

J45|, fl6 and antifcrromagnetic [16 
and ferromagnetic ordering |pl|,^4 
studies have been supplemented by the exact diagonal- 
ization method |f50| , [5l|| and by a number of approximate 
methods that are reliable in certain limits. Let us men- 
tion the non-crossing approximation (NCA) 
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as a strong-coupling approach, weak-coupling methods 
]5~I| , [Sp , |56| , f>l\ , variational approaches [^8[ |5{| |o) and 
others ]6l|, |62|, |63[| . Moreover, the d = oo model is phys- 
ically meaningful since the local approximation for the 
self-energy appears to be reasonable down to d = 3 Q 
or even d = 2 |35|. The essential physical properties of 
the d = oo Hubbard model are thus expected to be com- 
parable to those at d = 3 (d = 2). For the present study 
we therefore restrict ourselves to the d = oo case where 
we are able to judge of the reliability of a particular 
method on the firm basis that is provided by available 
QMC results. This should help us to come to conclusive 
results concerning the sum rule. 

Our main idea is to elucidate the meaning and the 
usefulness of the moment sum rule by considering a num- 
ber of standard approximations: (i) the simple Hubbard- 
I approximation (H-I) (ii) Hubbard's alloy-analogy 
solution (AA) |6q] , (iii) the so-called Edwards-Hertz ap- 
proach (EHA) [p7| in its improved version by Wermbter 
and Czycholl feq], and finally (iv) the generalization 
|69|| of the iterative perturbation theory |l5| to arbitrary 



band-fillings as proposed by Kajueter and Kotliar (KK) 
|7(i[) . It turns out that all approaches are inconsistent 
with the moment sum rule for m = 3 and yield the cor- 
rect moments only up to m = 2. 

The moments of the spectral density are intimately 
related to the coefficients in the high-energy expansion 
of the corresponding self-energy |7l| . We will show that 
thereby a possibility is opened to improve upon a partic- 
ular approximation analytically and to ensure that the 
moment sum rule is respected up to m = 3. 

This program has been completed successfully for 
the four approaches mentioned above. The result is an- 
other set of methods, each evolving straightforwardly 
from its above counterpart: (i) the spectral-density ap- 
proach (SDA) ]27|, || p2] |0), (ii) the modified alloy- 
analogy (MAA) [72|, 73|, (iii) the interpolating alloy- 
analogy-based approach (IAA) |24[], and (iv) the modi- 
fied perturbation theory (MPT) ]23|, 

Each of these new approaches can be considered as a 
conceptual improvement upon the original one: It can be 
shown that the additional inclusion of the fourth spec- 
tral moment does not affect their validity in different 
exactly solvable limiting cases. Furthermore, we will 
also show that significant improvement is achieved when 
comparing with the essentially exact QMC data of Refs. 

There is an interesting feature common to those four 
approaches that yield the correct moments up to m — 
3: The analytical expression for the self-energy depends 
on a higher-order correlation function B a which can be 
shown to be responsible for a (possibly) spin-dependent 
shift of the centers of gravity of the lower and the upper 
Hubbard band in the strong-correlation regime U i— ► oo. 
On the other hand, the original approaches are obtained 
if B a is (ad hoc) replaced by its atomic-limit (H-I, AA) 
or by its Hartree-Fock value (EHA, KK), respectively. 
One may thus expect that the inclusion of the fourth 
spectral moment is especially important what concerns 
magnetic order. 

After a general discussion of the moment sum rule 
and its conceptual importance in the next section, wc 
briefly discuss each of the methods in sections III- VI. 
Section VII presents new results for the paramagnetic 
phase on the hyper-cubic and for the ferromagnetic 
phase on an fec-type lattice. For the discussion a 
comparative opposition of the different methods against 
each other and against QMC results appears to be 
helpful. Our intention is to arrive at general conclusions 
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on the meaning and importance of the moment sum 
rule and of the correlation function B„ in particular 
which may be helpful for future analytical work. 



II. SPECTRAL MOMENTS AND 1/U 
PERTURBATION THEORY 

Let us start by recalling some essential facts for the 
limiting case of strong interaction (U 3> t). This will be 
important for the later discussion and for a concise for- 
mulation of the different approaches mentioned above. 
We first introduce the model as well as the basic physical 
quantities of interest. 

Using standard notations the Hubbard model reads: 



(1) 



The hopping integrals are assumed to be non-zero up 
to nearest-neighbor distances: Tuj\ = —t. The on-site 
energy Ta = To defines the energy zero. We consider an 
infinite-dimensional lattice with the usual scaling dt 2 = 
const, for the hopping to retain the model non-trivial 
||. Restricting ourselves to homogeneous phases (para- 
and fcrromagnetism) , the one-electron Green function, 



G a (E) = ((c ia ;cl)) E , 



(2) 



depends on the lattice geometry via the free Bloch den- 
sity of states (BDOS) only and can be written in the 
form: 



G a (E) 
where the BDOS 



dz 



(3) 



(4) 



is given in terms of the tight-binding dispersion 

e ( k )-^E^ k(R ' _R3) ^- (5) 

The expression ([|) for the Green function is based on 
the fact that the self-energy S cr (£') is k independent or 
site-diagonal for d = oo ^ [H?). Let us also introduce 
the k-resolved spectral density, 



A k<J (E) 



1. 



-Im- 



7T E + iQ+ - (e(k) - n) -Z a (E + i0+) ' 

(6) 



and the on-site or k-summed spectral density: 



The atomic limit of vanishing hopping (t = 0) rep- 
resents the zeroth-order result for the strong-coupling 
regime. We have: 



E a (E) = Un_ a 



U 2 n- a {l - Tl- a ) 
E + n - T - U(l - n- a ) ' 



(8) 



where n a — (f,\ a Ci a ). Let us now consider the first non- 
trivial correction to the atomic solution (||) in the strong- 
coupling regime. For U i— > oo perturbational results can 
be derived by performinga canonical transformation of 
the Hubbard model |^6[ [74| (75l (76) . As has been shown 
first by Harris and Lange [|26| , the 1/U perturbation 
theory predicts the spectrum to be dominated by two 
charge-excitation peaks (Hubbard bands); the weight 
of additional satellites in the spectrum is of the order 
(1/U) 4 and can thus be neglected for strong interaction. 
Furthermore, at each k point in the Brillouin zone the 
center of gravity T pa (k) as well as the spectral weight 
a pcr (k) of the lower (p = 1) and of the upper (p = 2) 
Hubbard band can be calculated exactly. Specializing 
the results of Harris and Lange to the d = oo case, one 
obtains: 

T la (k) = (l-n_ ff )e(k)+n_ ff B_ ff +OQ./U), 
T 2a (k) = (7 + n_ CT e(k) + (l-n_ CT )B_ CT + 0(l/[/), 
ai CT (k) = 1 - tl — (j 

+ - n_ CT )(B_ CT - e(k)) + 0{1/U) 2 , 



B„ 



To 



(9) 



r ^ I 3^yE^( c U.(2^-i)).(io) 



a 2(T (k) = 1 - ai CT (k) , 

where the B^ a is defined as 
1 



Contrary to the finite-dimensional case p^, Q , the 
k dependence of T pcr (k) and a p(T (k) is exclusively due to 
the Bloch dispersion e(k). For finite d an additional k 
dependence is introduced via B- a t— > Bk-<r- 

From the result (^|) valuable "global" information 
on the quasi-particle band structure can be read off: 
Apart from the Hubbard splitting, electron correlations 
manifest themselves in a specific narrowing of the lower 
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and upper Hubbard band, in a redistribution of spec- 
tral weight among them and in a shift of their centers 
of gravity that is given by n-„B- a and (1 — n-,j)B_ a , 
respectively. These (possibly spin-dependent) shifts, in 
particular, are expected to be important for ferromag- 
netic symmetry breaking: A spin-dependent shift of the 
band center of gravity may generate and stabilize ferro- 
magnetic solutions. 

Here the question arises how it can be ensured that 
a given analytical approach is consistent with the exact 
perturbational results of Harris and Lange for U i— > oo. 
For this purpose we consider the moments of the spectral 
density which for m = 0, 1, 2, . . . are defined as: 



M, 



(m) 
kcr 



E m A ka (E) dE. 



(11) 



that result from the above reasoning are the same as 
found by Harris and Lange in Eq. 

In most analytical approaches one obtains an (ap- 
proximate) expression for the self-energy. To make use 
of the sum rule, it becomes necessary to determine the 
spectral density from Eq. (||) and to calculate the mo- 
ments by carrying out the integration in Eq. (p]). Fi- 
nally, the (approximate) result for the moments has to 
be compared with the exact one given by Eq. (p"2|). 

The inconvenient integration in can be avoided 
when considering the high-energy expansion of the 
Green function and of the self-energy: In the spectral 
representation of the on-site Green function, 



(13) 



There is an alternative but equivalent way to calculate 
the moments: 



M, 



(m) _ 



(12) 



where CO — [O, if]_ denotes the commutator of an op- 
erator O with the Hamiltonian, and [•••,•••]+ is the 
anticommutator. Eqs. ([ll]) and ( |l2|) represent the mo- 
ment sum rule for the spectral density. In practice the 
sum rule is useful for low m only. The limitation arises 
from the fact that with increasing m Eq. (|l^) involves 
equal-time correlation functions of higher and higher or- 
der which are usually unknown. 

As has been shown in Ref. |33), the correctness of 
the moments up to m — 2 is a necessary condition for 
the evolution of the Hubbard bands as U i— ► oo. On 
the other hand, it is not sufficient as the following ex- 
ample shows: The self-consistent second-order perturba- 
tion theory in the interaction U for d = oo [ |54j predicts 
the correct moments up to m — 2 but does not yield the 
Hubbard splitting for large U |33|, Q . 

To be consistent with the results of 1/U perturbation 
theory, one has to check the existence of the Hubbard 
bands (which normally can be done easily). Then, it is 
sufficient to ensure that the moment sum rule is fulfilled 
for m — — 3: For each k point of the Brillouin zone, the 
first four moments provide four pieces of information on 
the spectral density that (for U > oo) unambiguously 
determine both, the centers of gravity and the weights, 
for the two Hubbard bands. Explicit expressions for the 
Hubbard model up to m = 3 can be found Refs. |2§| [35), 
for example. The dispersions T pa (k) and weights a pa (k) 



we expand the denominator in powers of 1/E. With Eq. 
(11) we immediately get: 



G4E) = hY, 



(m) 



rn—0 



(14) 



The coefficients in the high-energy expansion of the 
Green function are given by the (on-site) moments: 

Mi m) = M<™ } = N-^MgK Via Eq. f) they also 
determine the high-energy expansion coefficients of the 
self-energy: 



oo £,(m) 



m=0 



(15) 



Using the explicit expressions for the moments as ob- 
tained from Eq. (|l^) |28|, |35| , a straightforward calcula- 
tion yields: 



C ( 0) 
CP 



U 2 Tl- a (1 - TL — (j ) , 

U 2 n_ a (1 - tl — (j ) (B_ CT -n + U(l - n_ ff )) 



(16) 



Any analytical (approximate) expression for the self- 
energy can easily be checked against these rigorous re- 
sults simply by expanding in powers of 1/E. Provided 
that the coefficients ci™" 1 turn out to be correct up to 
m = 2, the moments of the resulting spectral density 
will be correct up to m — 3. For U i— ► oo this en- 
sures complete consistency with the 1/U perturbational 
results of Harris and Lange. 
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Finally, we have to express the expectation values in 
Eq. (|l^) in terms of known quantities. For the (spin- 
dependent) average occupation number we have: 



2 

h/3 



n=0 



(17) 



where E n — (2n+l)7r//3 and (3 = l/k^T. We also define 
the band- filling n — nj+nj. Ferromagnetic order is indi- 
cated by an asymmetry nj 7^ nj. in the spin-dependent 
occupation numbers. In this case the (dimensionless) 
spontaneous magnetization is given by m = m — ny. 
We can also determine B a (see Refs. |24|, p5| , for 

example) : 



1 2 00 / 2 
B a =T + — -— RcJ] ( -£ CT (i£„) - 1 



- n a ) hp 



[{iE n - ^ a {iE n ) - T + n) G a {iE n ) - h] 



(18) 



Fortunately it can be expressed in terms of the one- 
electron Green function and self-energy, despite the 
fact that the definition of B a ( |l0| ) includes higher-order 
correlation functions. 



III. TWO-POLE GREEN FUNCTION 

In the following we consider four different standard 
approximations for the self-energy: The Hubbard-I 
solution (Sec. Ill), the alloy analogy solution (Sec. IV), 
the Edwards-Hertz approach (Sec. V) and the ansatz of 
Kajueter and Kotliar (Sec. VI). Carrying out the high- 
energy expansion, we can check to which order m the 
moment sum rule is fulfilled. Furthermore, we will show 
up how the respective approaches can be corrected (if 
necessary) such that the sum rule is obeyed up to m = 3. 



A. Hubbard-I approximation (H-I) 

The so-called Hubbard-I solution is the simplest ap- 
proach that predicts a splitting of the non-interacting 
Bloch band into the two quasi-particle (Hubbard) sub- 
bands. In the original work Q it is derived by decoupling 
the hierarchy of equations of motion at the second level, 
i. e. by assuming {{ci a ni- a ;(X ja )) « (n. i _ CT )((c io -; c] CT )). 
From the resulting approximate Green function and 



from the Dyson equation one obtains the Hubbard-I self- 
energy which is identical with the self-energy (||) of the 
atomic limit. The high-energy expansion is readily per- 
formed, the expansion coefficients read: 



(7(0,H-I) 

£(1,11-1) 
r-(2,H-i) 



= U 2 n- 



(1 



C/ 2 n_ CT (1 - tl — (j ) (T -n+ U(l - tx — (j ) ) • 

(19) 



Comparing with (16), we notice that the m = 2 coeffi- 



cient (and thereby the m = 3 moment) is incorrect. 



B. Spectral-density approach (SDA) 

It is easily seen from the expansion coefficients that 
the sum rule up to order m = 3 can be restored if the 
atomic level T in the H-I self-energy (|8|) is replaced by 
T 1 ► B_ a . We obtain: 



S (SDA )( £ ) = [7n _ CT 



U 2 n- a (l - n-«r) 



E 



li- 



B- 



U(l 



(20) 

The Green function can be calculated from Dyson's 
equation, and from (18) we get B- a . We thus arrive 
at a conceptually simple improvement of the H-I solu- 
tion that introduces an effective atomic level which is 
determined self-consistently. 

It turns out that this is identical with the spectral- 
density approach (SDA) |^8|, |77], |35|, ||(| in infinite 
dimensions. The SDA is essentially equivalent to the 
Roth two-pole approximation for the Green function 
|7S| , [32]] . Furthermore, the SDA self-energy can also 
be obtained by means of the Mori-Zwanzig projection 
technique M, M, 



IV. ALLOY ANALOGY 

The main idea of an alloy-analogy solution for the 
Hubbard model is to consider the — a electrons to be 
"frozen" and to be randomly distributed over the sites 
of the lattice. The a electrons then move through a fic- 
titious two-component alloy that is characterized by two 
atomic levels E\ a and Ei a and the concentrations X\<j 
and X2a- The coherent potential approximation (CPA) 
can be used to perform the configurational average 
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over the positions of the frozen —a electrons. The self- 
energy for the a electrons is thus obtained via: 



o=5> 

P =i 



E pa — £,<j(E) — T 



pa 



l-lG a (E) [£ pCT -£ CT (£)-T ] 



(21) 



A. Hubbard's alloy-analogy solution (AA) 

Any two-component alloy analogy requires the spec- 
ification of the two atomic levels E pa and the concen- 
trations x pa . Within the conventional alloy-analogy so- 
lution due to Hubbard J6(| (AA) these are taken by re- 
ferring to the atomic limit: 



(AA) 



r (AA) 
X la 



E. 



(AA) 



To + U. 



.(AA) 



-a ■ (22) 

Using this ansatz in the general CPA equation ( pl| ) and 
rearranging the terms, we arrive at: 



S (AA) (jB) 



Un- 



l-lG a {E) (U-^ AA) (E)) 



(23) 



This AA solution turns out to be exact for a special but 
non-trivial limiting case of the Hubbard model: Switch- 
ing off the hopping of the —a electrons only (Ty i— > 
Tij a , Tij- a — SijTo), defines the Falicov-Kimball model 
(FKM) [^2| . As has been shown by Brandt and Mielsch 
pa] , in infinite dimensions the exact self-energy is given 
by Eq. @. 

To get the high-energy expansion coefficients of the 
AA self-energy, we insert the expansions (|l4|) and ( |T5| ) 
into the CPA equation (^lj) and sort the different terms 
in powers of 1/E. Considering all terms up to order 
l/E 2 yields the following set of equations: 



1 = 



= 



= 



= 



^ x pa ; 








^ Xpa (Ep a 

D 


-T -C<°>), 




{Epa 


-t -c(°>) 2 m(°)-c«' 


5Z X P& 


{Epa 


-T a -C^f(M^f 


P 

+ (Epa 


-T a 


_ C (0))2 M (1) 



2(E P<T -T - C^CWMf - C< a >l . (24) 



(0) 



Inserting the atomic levels and concentrations from (|2 
as well as the exact low-order moments M, 



1 and 



(i) 



Tq — fi + Uri-a and solving for C a results in 



^f(m,AA) (j( rn ) 



B- a 



*T 



_ £t(m,H-I) 



(25) 



for m = — 2. The high-energy expansion coefficients 
of the AA self-energy turn out to be identical to the 
coefficients within the H-I solution. Again the m = 2 
coefficient is found to be incorrect. 



B. Modified alloy-analogy (MAA) 

Simply replacing To i— > B- a in the CPA equation 
(21) and in ( p2| ) does not yield the correct m = 2 co- 
efficient since B—a would cancel out again in Eqs. (^4|). 
Another idea, however, turns out to be successful: The 
choice (|22] ) for the levels and concentrations is rather 
obvious. On the other hand, it is by no means predeter- 
mined. We therefore consider E ptT and x pa (p = 1,2) as 
free parameters to be fixed by the Eqs. ((2J) where ci™" 1 
and Mi-™'' are taken to be the exact coefficients (fj~6|). 
This automatically ensures the correctness of the mo- 
ments up to m = 3 and thereby provides an "optimized" 
alloy analogy. Solving Eqs. ( pi| ) for the parameters of 
the fictitious alloy yields: 



B (MAA) 
pa 



1 



[T + U + B^) + {-\y x 



i([/ + B_ CT -T ) 2 



Un_ a (T - B_ a ) 



and 



x 



(MAA) 



B- a + U(l - n_ CT ) - E 

(MAA) p(MAA) 

- 1(7 



(MAA) 



E. 



(26) 



(MAA) 
L 2a 



(27) 

It turns out that this result is identical to the result 
of the recently proposed modified alloy analogy (MAA) 
|72] , [73]| where it was derived by referring to the split- 
band regime @ of the CPA. 

Inserting the result into the CPA equation (|2|) yields 
the MAA self-energy: 



.(MAA) 



(E) = 



Un. 



j-Ga(E) (U-^ AA \E) 
1 - j-Ga(E)(B_ a -T ] 



(28) 



G 



One recognizes that the MAA reduces to the AA if B_ a 
is replaced by To. As can be seen from Eq. this is 
correct for the atomic limit and for the Falicov-Kimball 
model. Just as the AA, the MAA therefore remains 
exact within these two limits. It also reduces to the AA 
for a paramagnet at half-filling with symmetric BDOS 
where B-„ = Tq is required by particle-hole symmetry 



V. ALLOY-ANALOGY-BASED 
INTERPOLATION SCHEME 

A severe drawback of the alloy-analogy solutions con- 
sists in the fact that they are unable to reproduce the 
weak-coupling limit. This defect can be eliminated if the 
CPA equation is considered as a mere starting point for 
a reasonable interpolation formula which is demanded 
to be exact for small U as well as in the atomic limit 
and in the case of the FKM. 

Standard perturbation theory in the interaction U 
p5| provides us with the exact result for the self-energy 
in the weak-coupling regime. The first non-trivial per- 
turbational result for the self-energy, 

X a (E) = Un_ a + ^ oc \E), 



(29) 

beyond the Hartee term Un- a is given by the second- 
order contribution (SOC): 

^\x)A { ll(y)A ( ll(z) 



E (SOC) {E) 



IP 



An 



(f(x)f{-y)f(z) 



E — x + y — z 
f(- x )f(y)f{- z )) dx dy dz . 

(30) 

Here f(x) = l/(exp((3x) + 1) denotes the Fermi func- 
tion, and A a 1 \E) is the free (U = 0) spectral density 
being shifted in energy by a (possibly spin-dependent) 
constant E„: 



A^(E)=A^(E-E a ) 



(31) 



with A { ° ] {E) = -±ha.G a {E + i0+)\u= - The plain or 
conventional second-order perturbation theory (SOPT) 
is recovered for E a — 0, and with E a = U one 
obtains the SOPT around the Hartree-Fock solution 
(SOPT-HF) @ H, |1- To this end there is no need to 
specify the constant E a . The expression ([50]) is correct 
up to order U 2 for any function E a = E a (U) with 
S CT (0) = 0. 



A. Edwards-Hertz approach (EHA) 

An interpolation scheme that is based on the alloy- 
analogy idea but correctly accounts for the weak-coup- 
ling limit, has been suggested by Edwards and Hertz 
pTy . Within the Edwards-Hertz approximation (EHA) 
the self-energy has to be calculated from p7L p8[ : 



4 EHA )(£) = 



l-±G a (E) (l/-^ eha) (£)) 



(32) 



This differs from the conventional AA with respect to 
the propagator G a (E); The Green function G a {E) in 
Eq. ( p3| ) is (ad hoc) replaced by 



G a {E) 



U 2 n- 



1 



-Ei soc )(£-4 EHA )(£)+£ CT 



(33) 

where Ect SOC ' > (E) is the second-order contribution to the 
SOPT self-energy introduced in Eq. (|30|). A simple cal- 
culation shows that the EHA reproduces the atomic and 
the FKM limit for arbitrary band-fillings provided that 
the shifts E a in the definition ([u]) of the spectral density 
A^{E) and in Eq. ( |33| ) are determined self-consistently 
from the condition: 



f(E) A a (E) dE 



f(E)A£\E)dE. 



(34) 

This improvement upon the original theory |67| has been 
introduced by Wermbter and Czycholl |p8| . Expanding 
the EHA self-energy in powers of U up to order U 2 , one 
recovers the exact weak-coupling result (29). 

The high-energy expansion is performed straightfor- 
wardly. For the modified propagator one obtains: 



¥°^=e- + 



B 



(HF) 



/' 



r " 1 7^ 



(35) 



where Bi 11 ^ is the Hartree-Fock decoupled B a defined 
as: 

4 HF) =To+ ^"pT ^i^)' 1 ' (36) 
n a >{\-n)y>) jjti 



and (cLc lfT )W = -^JdEf(E) Im G ija {E + iO 



t c . \d) - _J_ 
E a ) I u=o- With this result at hand, the high-energy ex- 
pansion coefficients are obtained from Eq. (p2|): 



£-<(m,EHA) 



(37) 
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Again, the first two coefficients C a °^ and C^P are 
predicted correctly while the m = 2 coefficient turns 
out to be wrong. Only in the atomic and in the FKM 
limit as well as for the symmetric case n = 1 we have 
B_ a = B^=T . 



B. Interpolating alloy-analogy approach (IAA) 

The EHA interpolation formula can be improved 
such that the moments up to to — 3 are reproduced cor- 
rectly: The interpolating alloy-analogy approach (IAA) 
Q starts from the CPA equation (|2l]) with the Green 
function replaced by the modified propagator G a {E) 
which is given by Eq. (|33|): 



(IAA) 



E 



(IAA) 



l-\G a {E) [Eil AA) -E a (E)-T ] 



(38) 

The high-energy expansion of this equation is then used 
to determine the (unknown) atomic levels and concen- 
trations. The only difference with respect to the ex- 
pansion (|24|) of the general CPA equation is that the 
coefficients for m = 0, 1 have now to be inter- 

preted as the moments of the modified propagator. For 
m — 0, 1 they can be read off from Eq. (3^). Solving for 
the unknowns yields the following result: 



Et AA) =T a + Uu + B 



B 



(HF 



)+(-l) P x 



l( 



U + B_ a -B 



(hf; 



Un- a ( B^P -B^,. 



(IAA) 
L la 



B_„-B 



(HF) 



T + U(l-n_ a )-E\ 



(IAA) 



p(IAA) p(IAA) 
h 2a - h la 



1 — X 



(IAA) 
2a 



Inserting into (38) results in: 



££ AA >(£) 



Un- 



TjGa(E) ( 



u- 



^ AA) (E) 



(39) 



(40) 



1 - j-G a (E) (i?_ t 



B 



This differs from the MAA self-energy with respect 
to the modified propagator; also Tq is replaced by 



B 



(HF) 



in (| 



By construction the moments of the 



Furthermore, the theory is exact up to order U 2 and in 
the atomic and FKM limit. 



VI. ITERATIVE PERTURBATION THEORY 

While for the alloy-analogy-based theories the cor- 
rect weak-coupling behavior has to be enforced artifi- 
cially, it can be automatically accounted for within a 
diagrammatic approach. Furthermore, a diagrammatic 
ansatz appears to be attractive since it allows to con- 
struct thermodynamically consistent theories which are 
conserving in the sense of Kadanoff and Baym |88| and 
thus respect Luttinger's sum rule ^9|. The simplest con- 
serving approximation is given by Hartree-Fock theory. 
It recovers the to = and m = 1 moments only. One 
additional moment (to = 2) is correct within the self- 
consistent second-order perturbation theory (SC-SOPT) 
|54| , J56|. No improvement, however, is achieved by 
higher-order conserving approximations: On the con- 
trary, within the fluctuation-exchange approximation 
(FLEX) 1 90 and also for its analogue in d — oo [p5j, 
already the 1/E coefficient in the expansion of the self- 
energy and thus the m — 2 moment turns out to be 
incorrect. 

The (non-conserving) SOPT around the Hartree- 
Fock solution (SOPT-HF) @ || |7) yields the correct 
moments up to to = 1. For a paramagnet at half- filling, 
however, the moments are correct even up to to = 3 (in 
infinite dimensions). Furthermore, for t = SOPT-HF 
recovers the atomic-limit solution (||) and may thus pro- 
vide a reasonable interpolation between the weak- and 
the strong-coupling regime at n = 1 . 

In infinite dimensions the Hubbard model can be self- 
consistently mapped onto the single-impurity Anderson 
model (SIAM) || |§: 

#SIAM = £(£rf - V)ct c a + ^(tka ~ ^Wha^ka 
a ka 

+ U n a n- a + ^ v ka (cicika + al^cj) . (41) 

ka 

The on-site Green function and the local self-energy of 
the d = oo Hubbard model are identical with the im- 
purity Green function and self-energy provided that the 
conduction band dispersion eka and the hybridization 
strength Vka in the SIAM are chosen such that the self- 
consistency condition, 

l 



spectral density are correct up to to = 3 in the IAA. A a (E+fi) = E-{e. d -n)-H„{E)-h (G a (E)) 



(42) 
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for the hybridization function A a (E) — J2k^ka/(E ~~ 
eta) is fulfilled. 

Within the context of the symmetric single-impurity 
Anderson model (SIAM), SOPT-HF has been recog- 
nized to be quite well behaved Furthermore, as for 
the Hubbard model, SOPT-HF respects the moment 
sum rule up to to = 3 and the atomic limit of vanishing 
hybridization. The so-called iterative perturbation 
theory (IPT) [|l5) makes use of these advantages by 
combining SOPT-HF for the SIAM with the self-con- 
sistent mapping of the d = oo Hubbard model onto 
the SIAM. Indeed, IPT yields convincing results as 
has been proven by comparison with QMC studies 
p4| . The advantageous properties of the IPT 



43 



are lost, however, if one considers the non-symmetric 
case n^l. 

A. Kajueter-Kotliar approach (KK) 

Kajueter and Kotliar J70| proposed an interpolating 
ansatz for the self-energy of the SIAM which for arbi- 
trary fillings is correct in the weak-coupling and the 
atomic limit. For the symmetric case the Kajueter- 
Kotliar approach (KK) reduces to the usual IPT. The 
interpolative self-energy is given by: 



Y, a (E) = Un- 



(SOC) 



(E) 



(43) 



where n a — (c\c a ) is the average occupancy of the im- 
f sod 

purity site and {E) the second-order contribution 

within SOPT-HF for the SIAM. The coefficients a a and 
b a as well as a fictitious chemical potential Jl a that ap- 
pears in the definition of the Hartree-Fock spectral den- 
sity [cf. Eq. (£0|) 



— Im — 

7T E 



h 



A*<r - e d - A a (E + (j,)- Uri- 



(44) 

are considered as free parameters. fi a is fixed by requir- 
ing \i = /j\u=o + S(0), i. e. the Luttinger sum rule in 
d — 00 Q. For any choice of a„ and b„ the first two 
moments to = and to = 1 are correct. Within the 
KK approach, the parameters a a and b a are determined 
from the moment sum rule for m — 2 and the exact 
atomic-limit solution |7Q]. This yields: 



71-^(1 - n-g) 



(45) 



and 



U(l- 2n- a ) + \i a - fx 
U 2 n { 'l(l - n { H) 



(46) 



where rSP = Ti 1 



ff(E)A^(E)dE. 
This procedure, however, leads to an incorrect result 
for the m = 3 moment. Inserting ( |45| ) and (Efj) into the 
ansatz (|43"|), expanding in powers of l/E, and comparing 
with the exact result (|l^), we have 



m.KK) 



(47) 



up to to = 2. Here, within the context of the SIAM we 
have introduced the following definitions: 



ng (l- ng ) ^ Vto<0 ^ (2 "- 



1)> (48) 



and 



2n 



(i) 



V k M^) (1) (49) 



which are related to the corresponding definitions for 
the Hubbard model (|l(]) and ( |36| ) by means of the self- 
consistent mapping. As for the EHA the moment sum 
rule for to — 3 is violated since generally B^J^ 7^ B- a . 



B. Modified perturbation theory (MPT) 

It is possible to improve the Kajueter-Kotliar ap- 
proach such that the moments up to to = 3 are recovered 
correctly This constitutes the modified perturba- 

tion theory (MPT) |^5| . We start from the same ansatz 
(43) for the self-energy of the SIAM. Contrary to the 
KK approach, however, the parameters are fitted to the 
m = 2 and the to = 3 moment. Using 



.(SOC) 



C/ 2 n«(l-n«) 



U 2 n ( ll{\ 



E 



KB 



(HF) 



E 2 



(50) 



the high-energy expansion of the ansatz ( ff3[ ) is easily 
performed. Comparing with the exact result for the ex- 
pansion coefficients in Eq. dl9), we have to choose 



B_ r 



H - 



si H / ) +/I CT + C/(l-2n_ (T )) 



U 2 n ( ll{\ - n W a ) 



(51) 
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for the parameter b a to ensure the correctness of the 
moments. a a is still given by Eq. ([45|). 

We also modify the condition that fixes the param- 
eter Ji a . Analogous to the condition (|34|) used in the 
EH A and IAA, we demand: 



(i) 



(52) 



Choosing 'jl (T to enforce the Luttinger theorem implies 
the theory to be intrinsically limited to zero tempera- 
ture. The condition ( [52] ) is less problematic. It also 
allows to perform calculations at finite temperatures. 
Numerical results for T = have shown |^5| that a sig- 
nificant violation of the Luttinger sum rule occurs for 
fillings that considerably differ from half-filling only. 

The MPT reduces to the KK approach for the sym- 
metric case of a paramagnet at half-filling and sym- 
metric BDOS: Here we have B^J 1 ^ = B a = €d due 
to particle-hole symmetry. For small U the MPT re- 
covers the exact weak-coupling result up to order U 2 . 
By construction the moment sum rule is obeyed up to 
r7i = 3 which ensures the theory to be consistent with 
the strong-coupling perturbational results of Harris and 
Lange and also implies the validity of the MPT in the 
atomic limit for arbitrary filling. 

For the following discussion we refer to the KK 
approach as the original approach of Ref. |70|| but with 
the Luttinger theorem replaced by Eq. (|52j) as the 
condition to fix ~\x a . This allows to compare the "KK" 
approach with the MPT also for T ^ 0. 

VII. RESULTS AND DISCUSSION 

Numerical calculations have been performed for all 
approaches that have been discussed in the preceding 
sections. To study the effect of B a and the meaning of 
the moment sum rule for the paramagnet, we consider 
the hyper-cubic (he) lattice in infinite dimensions. The 



coordination number is Z = 2d. We take T, 



t 



t* j\J2Z (t* = const.) for the scaling of the hopping as 
d i— ► oo ||. The non- interacting Bloch density of states 
(BDOS) is a Gaussian §: 

1 



-(E/ey 



t* 



(53) 



To = Ta = has been chosen to fix the energy zero. 

We also consider a generalization of the d = 3 fee 
lattice to infinite dimensions 1 92 which favors ferromag- 
netic order 0,0. The hopping is scaled as t = t* j\/Z 



where Z = 2d(d - 1), and the BDOS reads: 

-(lW2£/n/2 

P (B \E)= 6 (54) 
Energy units are chosen such that t* = 1. 

A. Paramagnet 

Let us consider the he lattice first. In a wide region 
of the U -n plane the system is a paramagnetic Fermi liq- 
uid. Near half-filling antiferromagnetic order is observed 
H|, |l|, ||, [H HU ^ . Saturated ferromagnetism can 
be excluded completely |39| while non-saturated ferro- 
magnetic order has been found recently for very strong 
interaction U |53) . We will restrict our investigation for 
the he lattice to the paramagnetic phase only. 

According to the results of Harris and Lange, the 
first non-trivial effect of 1/U perturbation theory be- 
yond the atomic limit consists in a specific shift of the 
centers of gravity of the Hubbard bands. For the lower 
one it is given by n^ a B^ a . Fig. 1 shows this band 
shift as a function of the filling n for the he lattice 
at finite inverse temperature (3 = 7.2 and at moder- 
ate coupling strength U = 4 (the U dependence will 
be discussed below). For each of the approaches that 
obey the moment sum rule up to m = 3, the corre- 
lation function B^ a (inset in F ig. 1) has been deter- 
mined self-consistently via Eq. (18). Electron correla- 
tions are meaningless for filling n = 0; consequently, the 
band shift n- a B- a vanishes. At half-filling n = 1 we 
have n,- a B- a = due to particle-hole symmetry. In- 
between the band shift acquires positive values with a 
maximum at 77 « 0.5 — 0.75 depending on the approxi- 
mation scheme considered. The typical size of the shift 
(77_ (T i?_ cr 0.2) is significant compared with the vari- 
ance A = 1/V2 of the BDOS. 

For the same set of parameters but with the filling 
fixed at 77, = 0.79, Fig. 2 shows the quasi-particle den- 
sity of states (DOS) A a (E). For the discussion of the 
effects of the band shift we concentrate on the result for 
the spectral-density approach (SDA) and the Hubbard-I 
solution (H-I) first. Within the two-pole approaches the 
DOS for U i — > oo can be written in terms of the BDOS 



A a {E) 



,(B) 



E + n- 



jB—r, 
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p 



(B) 



E + fi-U-(l -n- a )B. 



,(55) 



where B- a i— > Tq for the H-I solution. Comparing the 
SDA with the H-I solution, the shift n^cB-c of the 
lower Hubbard band present in the SDA must be can- 
celed exactly by a corresponding shift of the chemical 
potential fj, for a fixed filling n. Thus the small differ- 
ences between SDA and H-I seen in the lower Hubbard 
band exclusively result from the finite value for U and 
will disappear for U i— > oo. On the other hand, the upper 
Hubbard band is significantly shifted to higher energies 
compared with the H-I solution. This effect persists also 
for U i— > oo where the energetic positions will differ by 



n B 




1 . < 




' 1 

MPT 


IAA 





0.10 



0.0 



0.5 

band filling n 



1.0 



o.oo 



Fig. 1. Filling dependence of the shift n_ CT B_ CT of the 
center of gravity of the lower Hubbard band as obtained 
self-consistently within the spectral-density approach (SDA), 
the modified alloy analogy (MA A), the interpolating alloy- 
analogy-based approach (IAA) and the modified perturba- 
tion theory (MPT). Inset: filling-dependence of B~ a . Lower 
panel: difference n- a {B- a - B^ F) ) for the IAA and MPT. 
Calculations for the d — oo paramagnetic Hubbard model 
on the he lattice at U = 4 and (5 = 7.2. The results are 
symmetric to the n — 1 axis (half-filling) . Energy units are 
chosen such that t* = 1. 



an amount (1 — 2n- a )B- a . 

Comparing the other pairs of approximation schemes 
among each other, we notice similar qualitative effects. 
Generally, the inclusion of B_ a in the theory results 
in fairly small changes of the lower but in a significant 
shift of the upper Hubbard band to higher energies (for 
fixed n). For the MAA there is an additional change in 
the peak heights compared with the AA. For the IAA 
and the MPT the shift of the respective upper Hubbard 
band against the EHA and the KK result turns out to 
be much smaller when compared with the MAA/ A A or 
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Fig. 2. Density of states A tT (E) as a function of energy for 
the he lattice at U = 4, n = 0.79 and = 7.2 as obtained 
within the Hubbard-I solution (H-I), the conventional alloy 
analogy (AA), the Kajueter-Kotliar approach (KK) and the 
Edwards-Hertz approach (EHA): dashed lines. Solid lines: 
resulting DOS within the improved theories: SDA, MAA, 
EHA and MPT. For comparison the numerically exact result 
obtained within the Quantum Monte Carlo method (QMC) 
by Jarrell and Pruschke pij is shown. 



11 



c: 

Z3 

i 

03 



LU 



1 1 ' 

SDA 


■ i 


| 1 i 1 i 


_MAA 






IAA 






~MDT 
Mr 1 






he 
" U=4 
- n=0.79 

P=7.2 

i 




i , i 



-2 2 4 6 
energy 

Fig. 3. Imaginary part of the self-energy as a function of 
energy for the he lattice at U = 4, n — 0.79 and j3 = 7.2. 
Thin horizontal lines indicate ImE(_E) = 0. Within the SDA 
ImE(£) consists of a single delta peak only which is indi- 
cated by a vertical solid line at the calculated energetic po- 
sition (also for the MAA). For the IAA ImE(_E) is always 
finite, large (negative) values near E = 2.2 are not shown. 



SDA/H-I. The reason is that within the EH A and the 
KK approach the band shift is not neglected completely 
(as for H-I and AA where n-uB-^ i— > n,- a To), but is 
taken into account on the Hartree-Fock level at least: 



. Indeed, as is shown in the second 



B^f^ ) as cal- 



pancl of Fig. 1, the difference rt- a (B- 
culated self-consistently within the IAA and the MPT is 
considerably smaller than the difference n_ CT (B_ (T — To) 
relevant for the SDA and MAA. 

Fig. 2 also shows the essentially exact result for 
the DOS as obtained by Jarrell and Pruschke H] by 
means of the QMC method. One clearly recognizes the 
two Hubbard bands which are well separated from each 



other. The energetic difference is approximately given 
by U. Additionally, in the vicinity of E = there is a 
rather sharp quasi-particle resonance which is reminis- 
cent of the Kondo peak in the SIAM. All approximate 
methods reliably reproduce the energetic positions as 
well as the weights of the Hubbard bands. If at all, 
there is only a slight improvement with respect to po- 
sitions and weights when taking into account the corre- 
lation function additionally. We can conclude that 
in this respect and in the case of the paramagnet the 
effects of B—fj are rather unimportant; considering the 
moment sum rule does not help much to improve the 
agreement with the QMC result. 

Since quasi-particle damping is neglected completely 
within the SDA, the Hubbard peaks turn out to be too 
narrow. Furthermore, there is no additional structure 
at E = 0. Obviously, the Kondo-type resonance cannot 
evolve in a two-pole approach. Somewhat broader peaks 
are predicted by the MAA; compared with the QMC 
density of states, however, their width is still too small. 
The IAA yields a considerable improvement. Yet, the 
size of the Hubbard gap is underestimated, and only 
a small asymmetry in the shape of the lower Hubbard 
band hints towards the resonance. The MPT achieves a 
reasonable agreement with the exact result. The Kondo- 
type resonance is clearly visible but not as pronounced 
as in the QMC spectrum. 

Additional insight is provided by the imaginary part 
of the retarded sclf-cncrgy which is shown in Fig. 3. The 
pole of the SDA self-energy © at E = B- a + U(l - 
ri—a) — fi gives rise to a single delta peak in Im £(£*), 
which, however, is meaningless since it falls into the 
Hubbard gap. Taking the imaginary part of the CPA 
equation ( |2l| ) at an energy E where ImG(E) = (is 
exponentially small), we obtain: 



= ImE CT (E)^ 



i 1 



\G a {E){E pa -T ) 



(56) 

This implies that a finite imaginary part of the MAA 
self-energy can be found in those energy ranges with 
non-vanishing DOS only, except for the energy point at 
which there is a zero of the denominator in Eq. (j56|). 
Here we find a delta peak in Im again (vertical 

line in Fig. 3). The delta peak does not contribute to 
the damping of the DOS but acquires most of the weight. 
This explains why the broadening of the Hubbard bands 
is still underestimated by the MAA. 
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In the IAA G a (E) is replaced by the modified propa- 
gator Go- (E) which yields a considerably stronger quasi- 
particle damping. We also note a small dip in the imag- 
inary part of the IAA self-energy at E 0. This 
tends to reproduce the correct Fermi-liquid behavior 
Im'E(E) ~ E 2 for T = 0. A T = calculation, however, 
shows ImS(O) ^ at U = 4 and n = 0.79. Similar 
to the EHA j67], ||], the IAA does not yield a Fermi 
surface in wide regions of the U-n plane. Contrary, the 
MPT self-energy always vanishes quadratically at E = 
and T = 0. Fig. 3 also shows a more pronounced dip 
in ImE(£) at E « in the finite-temperature MPT 
result. Compared with the IAA, the damping effect is 
of the same order of magnitude. In all cases the total 
weight / lmT, a (E) dE is given by — 7r[/ 2 n_ (J (l — 
This follows directly from the high-energy asymptotics 
of the self-energy (|l6|). 

There is a strong dependence of the band shift on 
the interaction strength U. Fig. 4 shows B_ a as a func- 
tion of U as resulting from the different approximation 
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Fig. 4. U dependence of the correlation function B- a (he 
lattice, n = 0.79, (3 = 7.2). For the IAA and the MPT 
-B^ F ' is shown additionally. 



schemes for n = 0.79 and /3 = 7.2. For large Z7 the cor- 
relation function B^ a becomes essentially independent 
of U. Within all approaches this asymptotic behavior 
is reached for U ~ 4; for U > 4 there is a very slight 
U dependence only. The asymptotic values predicted 
by the SDA and the MAA are considerably larger than 
those predicted by the IAA and MPT. Although the un- 
derlying physical concepts are quite different, there is a 
remarkable similarity between the results from the lat- 
ter two approaches over the whole U range. For small 
U the results are nearly identical (up to U w 2). This 
had to be expected since the IAA as well as the MPT 
self-energies are correct up to order U 2 . Via Eq. ( |lg| ) 
this implies that B^ a (U) is exact in order U, i. e. at 
U = B- a (U) starts with the correct slope within the 
IAA and MPT. Of course, at U = all approaches yield 
the same value for B- a . 

For the IAA and the MPT Fig. 4 also shows the de- 
pendence of B^J^ on U. Within the IAA B^f^ is given 
by Eq. @. The condition @ that fixes the shifts E a 
can be written as — n a . It implies that rig as 
well as (c^Cja)^ are independent of U, and therefore 
B_ a turns out to be a constant that depends of n and 
(3 only. Analogously, due to the condition (|5^) there is 
no U dependence of in the MPT. The hybridiza- 
tion function A CT , however, does depend on U due to 
the self-consistent mapping procedure. This leads to a 
U dependence of the second factor in Eq. (^) and thus 
of B_J' as can be seen in Fig. 4. 

For all approaches considered here we did not 
find ferromagnetic solutions on the he lattice. We 
calculated the uniform static susceptibility by applying 
an infinitesimally small external field: x = dm/dH. 
The susceptibility never diverged for U < 6. This is 
consistent with the QMC results of Refs. |l2|, and 
also with the complete instability of the Nagaoka state 
on the he lattice [^9). Let us mention, however, that a 
partially polarized ferromagnetic state has been found 
recently for extremely strong interaction U f53fl . 

B. Ferromagnetism 

While the moment sum rule for m = 3 has turned 
out to be of marginal importance for the paramagnetic 
phase, it is much more decisive with respect to the pos- 
sibility and the characteristics of ferromagnetic order. 
This shall be elucidated by the following arguments: 
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Firstly, ferromagnetism is a strong-coupling phe- 
nomenon. We thus have to account for the Hubbard 
splitting. It is a necessary condition for the evolution of 
the Hubbard bands as U > oo, however, that the mo- 
ment sum rule is obeyed up to to — 2 (see Ref. Q , for 
example). This is achieved e. g. by the H-I and the AA 
solution. 

Secondly, within the H-I approach ferromagnetism 
turns out to be rather unlikely. Ferromagnetic order 
may be possible for small n only when there is a large 
BDOS at the Fermi energy away from the band center 
of gravity [Q . The AA solution is even more prohibitive 
with respect to ferromagnetic order |)3|, Q . The reason 
|94| for this tendency to paramagnetism is that even 
in a (possible) ferromagnetic phase, the band centers 
of gravity T pcr (k) of the majority and of the minority 
electrons at each k point are equal: T p f (k) = T p y (k) for 
both, the lower (p = 1) and the upper (p — 2) Hubbard 
band. What is missing in the H-I and the AA approach 
to get a stable and extended ferromagnetic phase, is the 
possibility for a spin- dependent shift of the centers of 
gravity of the Hubbard bands. 

Thirdly, the possibility for a spin-dependent band 
shift is not only required but in fact is predicted by 
the 1/U perturbational approach of Harris and Lange 
which is exact for U i— > oo: As is seen from Eq. 
the spin-dependent shifts of the lower and the upper a 
band are n- a B- a and (1 — n^ a )B- a , respectively. The 
higher-order correlation functions included in the defi- 
nition @ of B— a open the channel for ferromagnetic 
solutions in the Hubbard model. 

Finally, as has been discussed in Sec. 2, the perturba- 
tional results of Harris and Lange for the strong-coupling 
regime can only be recovered if the moment sum rule is 
respected up to to = 3. 

It goes without saying that the 1/U perturbation 
theory merely determines the global features of the 
quasi-particle spectrum and abstracts from any internal 
structure of the Hubbard bands which may be due to 
damping effects or low-energy excitations, for example. 
Indeed, in the preceding section we have seen that differ- 
ent approximations yield rather different results even if 
all respect the sum rule up to to = 3. Thus it has to be 
expected that there is a dependence of the results on the 
approximation used also in the case of the ferromagnet. 

Our intention here is twofold: Firstly, as has been 
noticed in the case of the paramagnet, conceptually im- 
proving the method (SDA i-> MAA h-> IAA, MPT) also 
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Fig. 5. Filling dependence of the Curie temperature To for 
the SDA (solid line) and H-I (dashed line) as well as for the 
MAA (circles) in comparison with the QMC results (error 
bars) by Ulmke jl9| for the d = oo fec-type lattice. (Within 
the A A there is no ferromagnetic instability at all). 



yields a closer agreement with the QMC results. This 
will be checked for the ferromagnet, too. Secondly, com- 
paring the results of different methods with each other 
(H-I <-> SDA, AA <-► MAA, EHA <-► IAA, KK <-> MPT), 
we are able to "switch on" and to "switch off" in a con- 
trolled way the possibility for a spin-dependent band 
shift that is given via B_ a . This will clarify whether 
the concept of the spin-dependent band shift really helps 
to understand ferromagnetism, i. e. whether the system 
"uses" this possibility and "realizes" the spin-dependent 
shift of the bands. 

For the infinite-dimensional fec-type lattice with the 
BDOS given by Eq. ferromagnetic order even oc- 

curs at a moderate coupling strength U = 4 and extends 
over a rather wide region in the n-T plane. This was 
proven numerically within the QMC method by Ulmke 
Gj| . Due to the strong frustration of the fec-type lat- 
tice antiferromagnetic order is not expected. The highly 
asymmetric BDOS ( |5d| ) with its square-root singularity 
at the lower band edge E\, = — 1/V2 favors ferromag- 
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Fig. 6. Tc(n) as obtained within the IAA (solid line) and 
the EHA (dashed line) compared with the QMC result Q . 



Fig. 7. Tc(n) as obtained within the MPT (solid line) and 
KK (dashed line) compared with the QMC result H . 



netic order g|, g§ |9J. The QMC results of Ulmke @ 
are shown in Fig. 5. The filling-dependence of the Curie 
temperature Tq has been obtained from the zero of % 
which shows a linear Curie- Weiss behavior. 

The high density of states at the lower band edge is 
sufficient to produce stable ferromagnetism within the 
H-I solution. The resulting filling dependence of Tc is 
shown in Fig. 5. The ferromagnetic region in the n-T 
plane, however, is confined to very low densities only. 
As mentioned this is known to be typical for the H-I 
approach |2j. Contrary to the H-I solution, the SDA 
allows for a spin-dependent band shift. Fig. 5 shows 
that this leads to a considerable increase of the Curie 
temperature. Furthermore, ferromagnetic solutions are 
found for all fillings up to half-filling n — 1. These 
effects of -B-o- are much stronger than those found for 
the paramagnetic phase. 

At n = 0.58 and T = the system is fully polarized 
(m = n). The spin splitting of the center of gravity of 
the lower Hubbard band amounts to n^B^ — n^B^ = 
0.73. This has to be removed as T approaches the Curie 
temperature while the Hubbard gap still exists above 
Tc ■ The spin splitting is considerably smaller compared 



with the exchange splitting Un-t — Un^ = 2.32 within 
Hartree-Fock (Stoner) theory where correlation effects 
are neglected altogether. Consequently, the SDA yields 
a Curie temperature that is smaller by more than a fac- 
tor 3: /fc B lf DA) = 0.19 and k B T c Stoncr) = 0.74 (at 
n = 0.58, U = 4). The difference becomes larger and 
larger with increasing U since as a function of U the 
Curie temperature is unbounded in the Stoner theory 
while it reaches saturation within the SDA |3(J . The un- 
realistically high Stoner Curie temperature results from 
the necessity to bridge an exchange splitting of order U 
by the thermal energy k B Tc- 

A spin-dependent band shift is also realized within 
the MAA approach. Stable ferromagnetism is found for 
fillings < n < 1 (close to half filling we have not 
succeeded to obtain truely converged numerical solu- 
tions). Compared with the SDA, the Curie temperature 
is significantly lower for all fillings. While the mecha- 
nism leading to ferromagnetic symmetry breaking is the 
same in the SDA and MAA, ferromagnetism is destabi- 
lized to some extent by means of quasi-particle damping 



(see also Ref. 72 ). The broadening of the respective 
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spin-dependent bands tends to enhance their overlap 
and leads to a (self-consistent) depression of the mag- 
netization and thereby of Tc ■ The decisive meaning of 
the 771 == 3 sum rule becomes most apparent when we 
"switch off" the possibility of a spin-dependent band 
shift, i. e. if we compare the MAA results with those of 
the conventional AA. There is no ferromagnetic instabil- 
ity at all in the AA as has been proven in Ref. |)3| and 
also tested numerically for the fcc-type lattice considered 
here. Therefore, ignoring the possible spin-dependent 
band shift n- a B- a , misses the most important route to 
ferromagnetic order. 

Quasi-particle damping becomes stronger when turn- 
ing from the MAA to the IAA. Again, we notice a consid- 
erable decrease of the Curie temperature for all fillings. 
Fig. 6 shows that the maximum Tc in the IAA is less 
than half of the maximum Tq predicted by the MAA. 
Moreover, the ferromagnetic region shrinks to 77 < 0.61. 
Comparing with the EHA, we again notice that the 
spin-dependent band shift induced by B^ a favors fer- 
romagnetic order: Since is replaced by B^f^ in 
the EHA [cf. Eq. @)], the EHA partially accounts for 
the band shift. However, we find the spin splitting of 
B-cr to be larger than that of B^J^ . Consequently, Tc 
is enhanced in the IAA although the differences seen in 
Fig. 6 are moderate compared with the cases H-I/SDA 
and AA/MAA. 

The same arguments hold for the MPT and the KK 
approach: The spin splitting of the center of gravity 
due to B^J^ is somewhat smaller than that due to B- a . 
This implies a higher Tc in the MPT compared with the 
KK approach [see Eq. (Bt])]. The filling dependence of 
Tc is shown in Fig. 7. The differences between the KK 
approach and the MPT are moderate again and compa- 
rable to the EHA/IAA case. Once more, this suggests 
a one-to-one correspondence between the critical tem- 
perature and the band shift due to B_ a in the m = 3 
moment. 

For comparison we included the essentially exact 
QMC result || for T c {n) in Figs. 5-7. While the Stoner 
theory yields a Curie temperature that is by more than 
one order of magnitude too high (see Ref. [p5[ , for ex- 
ample) , realistic values are found by the approximations 
discussed here. Even the simplest scheme consistent 
with the requirements of Harris and Lange, the SDA, 
considerably improves upon the Stoner result. A quan- 
titative agreement with the essentially exact result, how- 
ever, cannot be expected since damping effects are ne- 



glected completely. A further improvement towards the 
QMC result is achieved by the MAA which takes into 
account a finite quasi-particle life time. Yet, the Curie 
temperatures are systematically too high since damp- 
ing effects are underestimated (cf. previous section and 
Fig. 9). The correct order of magnitude for Tc can only 
be expected if (i) the theory respects the m = 3 sum rule 
and (ii) realistically accounts for quasi-particle damping 
as is done by the IAA and the MPT. 

Fig. 8 focuses on the temperature dependence of the 
magnetization at 77 = 0.58. The QMC data can be well 
fitted to an S — 1/2 Brillouin-function. Extrapolation 
to T — yields m — n p9| ]. A fully polarized ground 
state is also obtained by the SDA and MAA, the MPT 
predicts a slightly lower 777. The filling 77 = 0.58 is near 
to the quantum-critical point ni IAA ' ) = 0.61 in the IAA 
(cf. Fig. 6). Since phase transitions have been found to 
be always of second order within the IAA, a partially 
polarized ground state near n£ IAA ^ is plausible from the 
result in Fig. 6. The vicinity to 77c IAA ' = 0.61 may 
also explain the non-Brillouin-function-like shape of the 
magnetization curve. Fig. 8 shows second-order transi- 
tions for the SDA, MAA and MPT at 77 = 0.58. De- 
pending on the filling, however, all three methods may 
also predict first-order transitions to the paramagnetic 
phase. Non-continuous transitions are observed for fill- 
ings larger than ~ 77 max where n max is the filling at which 
Tc reaches its maximum. For the MPT the change of 
second-order into first-order transitions is marked by the 
kink in the Tc(n) curve at 77 = 0.67 (Fig. 7). While mag- 
netic first-order transitions may seem to be implausible, 
they cannot be ruled out by a rigorous argument. On 
the other hand, we cannot strictly exclude that they are 
artefacts of the approximations. 

For the second-order transitions shown in Fig. 8 the 
homogeneous, static susceptibility \ = dm/dH\H=o 
diverges at T = Tc- Within all approaches consid- 
ered x obeys the Curie- Weiss law for high tempera- 
tures. But also for temperatures close to Tc we ob- 
serve an almost linear trend of x • This is consistent 
with the QMC result. Note that the different slopes in 
the plot result from the different Curie temperatures. 
From x — C / (k^T — fee©) we have calculated the (di- 
mensionless) Curie constant C. For the paramagnetic 
Curie temperature we have assumed O = Tc except 
for the IAA where fee© = 0.057 has been obtained by 
extrapolation of the linear trend of x~ 1 a t high temper- 
atures. The QMC value C = 0.47 is to be compared 
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Fig. 8. Homogeneous magnetization m as a function of the 
reduced temperature T/Tc for the fcc-type lattice at U — 4 
and n — 0.58 as obtained self-consistently within the SDA, 
MAA, IAA and the MPT (solid lines below T/T c = 1.0). 
Bottom panel: QMC data from Ulmke [0] (diamonds) . The 
dashed line below T/Tc = 1.0 is a Brillouin-function fit to 
the data. For T/Tc > 1-0 the inverse homogeneous static 
susceptibility x -1 (multiplied by a factor 4) is shown. Cir- 
cles: QMC data |49|. Dashed line: linear fit to the data. 
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with C = 0.42,0.52,0.50,0.57 for the SDA, MAA, IAA 
and MPT, respectively. 

Finally, we compare the spin-dependent spectral den- 
sity as obtained by the approximate schemes with the 
essentially exact QMC result. Fig. 9 shows A a (E) in 
the ferromagnetic phase for a =\ (majority) and a =\ 
(minority spin) at U = 4 and n — 0.58. For a mean- 
ingful comparison the respective temperature has been 
chosen to yield a magnetization of m = 0.4. In the 
SDA spectrum we recognize the large Hubbard split- 
ting ~ U and a much smaller additional spin splitting 
for each of the two Hubbard bands. Furthermore, a 
correlation-induced spin-dependent narrowing as well as 
a spin-dependent weight of the subbands is observed. 
All features agree qualitatively well with the predictions 
of Harris and Lange which are valid for U i— > oo. A 
more quantitative estimate shows up discrepancies with 
respect to the positions and weights given by Eq. (^) 
which are due to the moderate value assumed for U. 

The shape of the peaks only slightly deviates from 
the shape of the BDOS. This is an artefact resulting 
from the complete disregard of quasi-particle damping. 
Turning to the MAA spectrum, a broadening of the 
peaks is introduced. Damping effects are even more 
pronounced by the IAA and the MPT. This turns out 
to be quite realistic when comparing with the QMC 
spectrum which also shows up a considerable peak 
broadening. In all approximations the spin-dependent 
position and weight of the lower Hubbard band is well 
reproduced, the position of the upper Hubbard band, 
however, is shifted to higher energies compared with 
the QMC spectrum. This indicates that for U = 4 the 
interaction is not sufficiently strong to reproduce the 
position resulting from the l/U perturbation theory. 
An additional peak in the a =| as well as in the a =[ 
channel of the QMC spectrum shows up at E ss 0.8 the 
physical origin of which is unclear. It is not reproduced 
by the approximate methods. At E w there is also 
a weak shoulder in the a =\ QMC spectrum. This is 
clearly reproduced by the MPT; but also in the IAA 
spectrum there is a (very weak) structure at E m 0. 
As for the paramagnet the peak is interpreted as a 
Kondo-like resonance. It cannot show up in the J, 
spectrum, since the filling of the j band is too small. 

VIII. CONCLUSIONS 

The present paper has discussed the moment sum 



rule as a valuable source of a priori information on 
the physical properties of strongly interacting lattice 
fermion models. The analysis has been restricted to the 
single-band Hubbard model in infinite dimensions since 
this allows to compare with numerically exact solutions. 
Here the reliability of approximate approaches as well 
as the usefulness of the sum rule for practical improve- 
ments can be estimated safely. The general trends that 
show up for d — oo are expected to hold also for finite di- 
mensions, where extensive use of the sum rule has been 
made in the past. Let us briefly list up the conclusions 
evolving from our analysis: 

(i) The moment sum rule fixes the norm, the cen- 
ter of gravity, the variance, etc. of the interacting spec- 
tral density at each k point in the Brillouin zone. To 
some extent this determines the "global" single-particle 
excitation spectrum for arbitrary U. The sum rule is 
of special conceptual importance in the strong-coupling 
regime. We could argue that any approximation that 
predicts the Hubbard splitting for U i— ► oo and that 
yields the correct moments up to m = 3 is fully consis- 
tent with the first non-trivial results of the perturbation 
theory in l/U by Harris and Lange. This includes im- 
portant correlation effects such as the U, filling and spin 
dependence of the centers of gravity, the widths and the 
weights of the two Hubbard bands. 

(ii) The high-energy expansion of the Green func- 
tion and of the self-energy provides a practicable way 
to check to which order m a particular approximation is 
consistent with the sum rule. The majority of analyt- 
ical but approximate approaches commonly used turn 
out to be at variance with the sum rule for m = 3. We 
considered four methods in detail: The Hubbard-I and 
the Hubbard-III alloy-analogy solutions as well as the 
Edwards- Hertz and the Kajueter-Kotliar approaches. In 
each case the sum rule is fulfilled up to m — 2 only. 

(iii) We have shown that the compatibility with the 
m = 3 moment sum rule can be restored by slight mod- 
ifications of the original methods. Starting from the H-I 
and the conventional AA solution, this yields an opti- 
mized two-pole and an optimized alloy-analogy approach 
which turn out to be identical with the spectral-density 
approach and the so-called modified alloy-analogy solu- 
tion, respectively. Analogously, the interpolating alloy- 
analogy-based approximation and the modified pertur- 
bation theory evolve straightforwardly from the EHA 
and the KK approach. We have checked that the ap- 
plied modifications do not affect the validity of the orig- 
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inal methods in exactly solvable limiting cases of the 
Hubbard model. 

(iv) All approximation schemes have been imple- 
mented for numerical evaluation. We have presented 
and discussed new results for correlation functions and 
dynamical quantities which are directly compared with 
essentially exact QMC results on the hyper-cubic and 
an fcc-type lattice in infinite dimensions. We observed 
that with increasing complexity (which is necessary to 
recover an increasing number of exact limits) the agree- 
ment with the QMC results is improved. This corrob- 
orates the usefulness of an interpolative analytical ap- 
proach to the Hubbard model. 

In particular, we found it necessary to account for 
damping effects in a realistic way: If most of the total 
weight of the imaginary part of the self-energy does not 
contribute to quasi-particlc damping, as in the MAA 
and more extremely in the SDA, the peak widths that 
are obtained by QMC cannot be recovered. Compar- 
ing the results of the different methods we also found 
that quasi-particle damping has a destabilizing effect on 
ferromagnetism which manifests itself in an overall de- 
crease of the Curie temperature. 

(v) Out of all approximation schemes considered 
here, the modified perturbation theory turns out to be 
most reliable. In the paramagnetic phase it correctly re- 
covers the Kondo-type resonance, shows up the expected 
Fermi-liquid properties and approximately fulfills the 
Luttinger sum rule. The MPT yields the qualitatively 
correct density of states in the weak- and intermediate- 
coupling regime as confirmed by comparison with the 
QMC results. At half filling the approach reduces to 
the iterative perturbation theory (IPT) which is known 
to yield a rather realistic description of the Mott tran- 
sition. For the ferromagnetic phase the MPT is able 
to predict the Curie temperature within an accuracy of 
~ 30%. Thereby, the MPT turns out to be superior com- 
pared with the IAA. The main drawback of the latter 
surely consists in its inability to predict a Fermi sur- 
face in wide regions of the phase diagram. The same 
defect is found in the MAA. However, the straightfor- 
ward and physically well motivated concept of an "opti- 
mized alloy-analogy" remains rather attractive. Finally, 
the SDA is surely too simple to compete with the more 
elaborate methods. On the other hand, its simplicity is 
advantageous when a quick though rough estimation of 
the magnetic phase diagram is required. 

(vi) The comparative opposition of the different 



methods has particularly shown the importance of the 
m = 3 moment sum rule with respect to ferromag- 
netism. While the differences between the methods 
with and without regard of B a are unimportant or even 
negligibly small for the paramagnetic phase, the cor- 
relation function B a considerably affects the ferromag- 
netic/paramagnetic phase boundary as well as the Curie 
temperature. The effect is most striking when compar- 
ing the AA, where ferromagnetic instabilities are ruled 
out completely, with the MAA which yields a reasonable 
dependence Tq (n) for the fcc-type lattice and even over- 
estimates Tc- The reason is found in the term B a that 
appears first in the expression for the m = 3 moment. 
It opens the possibility for a spin- dependent shift of the 
centers of gravity of the Hubbard bands. In the AA 
and also in the H-I solution this is missing completely 
(B a i — > To) • In the EH A and the KK approach it is 
only approximately accounted for (B a t— > B^ HF ^). Since 
£>ct HF ^ is found to be too small in the EHA and the KK 
approach, both methods yield a lowered Tc- 

From the recent QMC studies plj, p[ in par- 
ticular (rum P. of. [!)")]. but also from variational resulls 
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the importance of the lattice structure for 
ferromagnetism well established. It is striking that all 
presently considered approximations that account for 
B a yield a ferromagnetic instability for the non-bipartite 
fcc-type lattice while (at least up to U < 6) the param- 
agnet is stable on the he lattice. In Ref. |Q Vollhardt 
et al. argue that a strongly asymmetric BDOS with high 
spectral weight at the lower band edge favors ferromag- 
netism because this minimizes the kinetic energy of the 
fully polarized state. The shape of the non-interacting 
BDOS is relevant since the f DOS is unrenormalized if 
n-f = n < 1. If it is assumed that details of the BDOS arc 
not so relevant for the paramagnetic state where corre- 
lation effects dominate, the ferromagnetic state becomes 
energetically more favorable. Considering the moment 
sum rule or, equivalently, the Harris and Lange results, 
the argumentation can be refined: For the fully polar- 
ized state the expression ( |Io|) for the band shift reduces 



to (1 — n)nSf 



--©kin- This implies 



that a strong asymmetry of the BDOS not only mini- 
mizes the kinetic energy (at a given filling n) but also 
maximizes the shift n^B^ > of the center of gravity 
of the lower { Hubbard band to higher energies. There- 
with, the spin splitting n^B^ — nj,Bj, = rz.j-.B-f of the 
lower Hubbard band becomes maximal which enhances 
the stability (reduces the instability) of the fully polar- 
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ized state. This mechanism is already present within the 
simple SDA. While the approach overestimates the pos- 
sibility for ferromagnetism, it predicts the correct trend 
when comparing different lattices |22| . 

As an intermediate to strong-coupling phenomenon, 
ferromagnetism in all its different aspects is hardly 
amenable to a simple explanation. We feel, how- 
ever, that by considering the spin-dependent band 
shift appearing in l/U perturbation theory and the 
m = 3 moment one is able to capture some of the 
essentials in a rather simple and physically intuitive way. 
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